Evidence of early multi-strange hadron freeze-out in high energy nuclear collisions 
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Recently reported transverse momentum distributions 
of strange hadrons produced in Pb(158AGeV) on Pb colli- 
sions and corresponding results from the relativistic quantum 
molecular dynamics (RQMD) approach are examined. We 
argue that the experimental observations favor a scenario in 
which multi-strange hadrons are formed and decouple from 
the system rather early at large energy densities (around 1 
GeV/fm 3 ). The systematics of the strange and non-strange 
particle spectra indicate that the observed transverse flow de- 
velops mainly in the late hadronic stages of these reactions. 



The purpose of the current and forthcoming heavy 
ion programs at the high-energy laboratories CERN 
(Switzerland) and Brookhaven National Laboratory 
(USA) is to probe strongly interacting matter under ex- 
treme conditions, i.e. at high densities and temperatures. 
The central subject of these studies is the transition from 
the quark-gluon plasma to hadronic matter. In the early 
phases of ultrarelativistic heavy ion collisions, when a 
hot, dense region is formed in the center of the reac- 
tion, there is copious production of up, down, and strange 
quarks. Transverse expansion is driven by the numerous 
scatterings among the incoming and produced particles. 
As the medium expands and cools, the quarks combine 
to form the hadrons that are eventually observed. 

In this Letter we are going to address the question 
whether data from heavy ion experiments allow one to 
make statements about the existence of strange hadrons 
in a medium of high energy density e > lGeV/fm . So 
far, experimental information about survival or "melt- 
ing" of quark bound states was restricted to charmonium 
states only |IJ. On the other side, the strange quark mass 
is intermediate between charm and the very light flavors 
(up and down). Strangeness should in fact not be con- 
sidered as an "external probe" of the hot medium like 
charm, as it may influence fundamental characteristics of 
the QCD phase transition itself, e.g. its order Fur- 
thermore, we argue in this Letter that the strange and 
non-strange hadron spectra encode information on the 
timing of flows and pressures. It has been recognized for 
a long time || that hadron momentum spectra are a valu- 
able source of information on the collective flow develop- 
ing in ultrarelativistic heavy ion collisions. The flows 
may be related to bulk and transport properties in the 
ultra-dense matter like transient pressure and viscosities. 
Presence of hydrodynamic behavior is expected at least 



for truly large colliding system such as Pb+Pb. Con- 
cerning the search for the quark-gluon plasma one would 
like to identify the regions in energy density at which 
the Equation of State (EOS) softens, presumably due to 
the phase transition or a cross-over between hadronic and 
quark matter. Another topic of interest is to see the EOS 
becoming hard again at yet higher densities. Asymptot- 
ically, the EOS approaches the Stefan-Boltzmann limit 
p = e/3, where p is the pressure, according to recent lat- 
tice calculations Q. Note, however, that the regime of 
perturbative quarks and gluons seems to be reached only 
at grand unification scales In the context of ultra- 
relativistic collisions the EOS dependence on energy and 
baryon density translates into a pressure dependence on 
time, because expansion dilutes the matter continuously. 
Here hadrons with varying strangeness content play an 
important role. With vastly differing reaction rates in the 
medium, they decouple at different times from the evolv- 
ing system. Pictorially speaking, we may employ these 
spectra to get a sequence of snapshots of the transverse 
flow present at each of the species-dependent decoupling 
times. 

On the experimental side, the predicted presence of 
strong radial transverse flow in the Pb(158AGeV) on Pb 
collisions (0.4 to 0.6 c) || has been deduced from the 
systematics of non-strange particle spectra already some 
time ago |7|J^]. The long awaited spectra of multiple 
strange hadrons <E>, S, and f2, measured at mid-rapidity, 
were reported during the Quark Matter '97 conference 
|j],[lO| . It came somewhat as a surprise that the reported 
slopes of these multi-strange hadrons (EE's and fi's and 
possibly <i>'s) are much softer than expected from the 
trends in the previously observed mass dependence of 
the slope parameters [0 . 

In this Letter we will argue that these deviations find 
their natural explanation in a reaction scenario in which 
the multi-strange hadrons freeze out early, before sizable 
transverse flow has been developed. The corollary state- 
ment is that the transverse flow is a phenomenon emerg- 
ing only later in the Pb on Pb reactions. In order to go 
beyond a qualitative interpretation we utilize a transport 
theoretical approach - relativistic quantum molecular dy- 
namics (RQMD) - whose predictions agree well with the 
observations. 

Several scenarios for the unknown earlier stages in ul- 
trarelativistic nucleus-nucleus are imaginable and have 
been put forward. We will try to use the most recent in- 
formation on the spectra to eliminate some of them. Let 
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us list some of the possible choices: 

• Is the expanding matter made of weakly interacting 
quarks and gluons, or of bound states (hadrons)? 

• Does the flow appear already at an early stage of 
the collision |ff|] or rather later Jl^]? 

• Does freeze-out occur instantaneously jl3| or se- 
quentially @fL§? 

We do not believe that the existence of strong trans- 
verse flow in the Pb(158AGeV) on Pb reactions belongs 
to the unresolved questions. It is well-known that the 
difference between flow - a space-momentum correlation 
- and excitations purely in momentum space (tempera- 
ture or random kicks |lq| ) shows up in observables which 
are sensitive to the phase-space densities of the finally 
emitted hadrons (like nucleon cluster formation |(| and 
HBT lUl). 
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FIG. 1. Measured slope parameters as a function of par- 
ticle mass. Preliminary results of the slope parameters of 
strange particles are shown with open symbols. The line rep- 
resents the parameterization mentioned in the text. 



Let us turn now to the recent data and interpret them 
in the light of these ideas. A compilation of the exper- 
imental data is shown in Fig. I. The non-strange data 
may be parameterized by T=Tf D + rn(/3 t } 2 , where T the 
slope parameter at mid-rapidity [jf8f and m the particle 
mass. The parameter values Tf «145 MeV, (flt) = 0.4c. 
Recent reports show that this trend is not only obeyed by 
7r + /~ ,p,p and deuterons but also by hadrons containing 
one single strange quark p 10 1. However, note the strong 
deviation of the multiple strange baryon (H - , Q) slopes 
from the general trends. How much the $ deviates from 
the general systematics is not yet settled [[f9| . 

We may draw three conclusions from the compilation 
of slope parameters. The frequently employed picture - 
one-fluid flow until break-up of the matter at a common 
freeze-out state - is untenable in view of these data. Such 



a picture results in a unique dependence of slope parame- 
ters on the particle mass, in clear contrast to the observa- 
tions. Note that feed-down corrections which could spoil 
the universal dependence are of relevance mostly for pi- 
ons but not the "heavy" hadron sector (with masses close 
to and above 1 GeV). Furthermore, the striking differ- 
ence between multi-strange hadron and the other hadron 
slopes points to dynamics in the hadronic stages that 
causes the non-universal pattern. It is well-known from 
the early ideas of "strangeness as a QGP signal" |20| that 
strange flavor impacts reactions in quark matter rather 
differently than in a hadron gas. For instance, perturba- 
tive QGP interactions at relevant temperatures are only 
mildly affected by the light quark mass variations (5- 
f60 MeV). On the other side, transport properties of the 
heavy baryons with its mass of 1672 MeV are expected 
to differ completely from the almost massless pions (140 
MeV) pl| . Finally, we note that the 5 slope parameters 
in Pb on Pb collisions are comparable to those in p on 
W and S on W reactions p2] ]. This fact provides strong 
constraints on any reaction dynamics. While much larger 
energy densities are surely produced in the heavy-ion col- 
lisions it does not seem to affect the transverse motion 
of particles like the H. Scenarios in which the precursor 
matter experiences sizable flow before freeze-out of the 
multi-strange hadrons seem therefore to be ruled out by 
the data j23j. The conclusion that there can be only a 
tiny collective component in the transverse momentum 
spectra of the multi-strange hadrons is corroborated by 
an analysis of the chemical composition. The slope pa- 
rameter 237 ± 24 MeV of the fl particle is consistent with 
the "temperature" value extracted from the particle ra- 
tios [El. 




FIG. 2. 
parameters. 
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The RQMD 2.3 model prediction of the slope 



We now turn to the RQMD calculations in order to go 
beyond qualitative statements. The RQMD model |25| ] 
provides a microscopic description of heavy ion collisions 
which has been highly successful in predicting most of 
the observed features over a wide range of conditions. We 
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have generated 1600 events for Pb(158AGeV) on Pb em- 
ploying RQMD (version 2.3) and analyzed the final spec- 
tra in the same fashion as has been done with the mea- 
surements. Before we analyze the aspects of the RQMD 
evolution dynamics which are pertinent for our discus- 
sion we present the results for the slope parameters of 
the various species in Fig. 2. Very good agreement is 
found between RQMD predictions and preliminary data. 



RQMD(v2.3) 158AGeV Pb+Pb (b<1 fm) 
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FIG. 3. Time evolution of the transverse source size for 
mid-rapidity baryons and pions. 



The transverse expansion dynamics generated by 
RQMD may be schematically decomposed into two 
stages. The pre-equilibrium stage is determined by the 
initial excitation and fragmentation of color strings and 
ropes followed by kinetic equilibration in an ultradense 
hadron gas. During this stage the effective transverse 
pressure is ultrasoft which is a result of the combined 
effects from initial collision geometry string dynamics 
which keeps a memory of the original beam axis, finite 
mean free paths and creation of "resonance matter" jl2) . 
According to the RQMD calculations this stage lasts for 
about 5 fm/c in central Pb on Pb reactions. The state 
when local kinetic equilibrium is finally achieved is soon 
followed by a break-down of equilibrium due to the di- 
luteness of the hadron gas and the finite size of the sys- 
tem. It has also been shown that the RQMD evolution 
of the multi-component hadronic fluid is characterized by 
non-ideal effects, even in the dense regime |lj|. The pions 
accelerate quite easily, their motion more or less governed 
by their own EOS. In contrast, the rarer heavy particles 
cannot keep up with the pion "fluid" and are left behind. 
The developing flow of matter according to the RQMD 
calculations is illustrated by plotting the mean transverse 
distance of hadrons from the center of the collision region 
(see Fig. 3). We see from Fig. 3 that the baryon matter 
does not expand at all during the first 5 fm/c. Only after 
the soft stage has elapsed the baryons develop collective 
flow which is reflected in the increase of their spatial dis- 
tribution. Essentially all of the transverse baryon flow 



is created during the late stage. Owing to the differ- 
ent transverse sizes of the collision region this stage lasts 
much longer in Pb on Pb than in reactions with smaller 
projectiles and leads to the marked difference between 
momentum spectra for Pb+Pb as compared to p+A and 
S+A reactions. 



RQMD(v2.3 cd) 158AGeV Pb+Pb (b<3fm) 
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FIG. 4. The time and transverse radius distributions of 
midrapidity Q's and nucleons at freeze-out in central Pb+Pb 
collisions at 158AGeV from RQMD. 

Why are the multiple strange particles not dragged 
with the heavy particle flow? Their interactions in the 
expansion stage are dominated by resonance formation. 
A good measure of the reaction rates are therefore the de- 
cay widths of the baryon resonances to which the baryons 
couple. We see from the particle data group tables that 
the decay widths of resonances is a strong function of 
their flavor content ||(J. Approximately, the trend is 
0.45:0.62:0.85:1 for fl*:E*:Y*:N* . We expect from these 
numbers that the 5 collision rates will be suppressed by 
on the order of 30 to 40 percent compared to A's and 
nucleons. Furthermore, the is basically not involved in 
forming these resonances. A reason is that the ir-il sys- 
tem does not match any of their flavor quantum numbers. 
Furthermore, fl+ij and collisions are suppressed 

due to the comparably rarer meson partners. The physics 
of resonance formation is built into the RQMD approach 
(for details see p5| ). RQMD predicts considerably ear- 
lier freeze-out of Q baryons (and S's) than nucleons (and 
A's). The time and spatial distributions of nucleons and 
fi's are displayed in Fig. 4. Most of the fi's freeze-out 
between 2 and 8 fm/c while the nucleon freeze-out distri- 
bution is centered around 20 fm/c. Note that the freeze- 
out time spectra are characterized by sizable differences 
between median and averages. Very soft collisions be- 
tween particles of almost equal velocities lead to the long 
tails in the distributions but are of no importance for the 
reaction kinetics. 

The corresponding transverse distance distribution 
which is also shown in Fig. 4 reveals that the f2 source at 
freeze-out is very similar to the initial source. In contrast, 
nucleons are transported by the collective flow into outer 
shells. We can analyze the RQMD results concerning the 
local energy densities at which the multi-strange hadrons 
decouple from the system. Typically, these local energy 
densities cluster around 1 GeV/fm 3 . Clearly, this sets a 
lower limit for the densities at which these particles may 
form as quark bound states. 
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We expect that the RQMD evolution in the late dilute 
hadron gas stage is calculated rather accurately - under 
the assumed given initial conditions. For this stage (and 
only here) our inability to calculate nucleus- nucleus reac- 
tions based on quantum chromodynamics (QCD) can be 
compensated. Utilizing the empirical information about 
the interactions between hadrons, kinetic equations can 
be set up and solved like in RQMD and other approaches 
p5fl , p7-p9|]. Theoretical justification for semi-classical 
transport comes essentially from the particle-like behav- 
ior with the DeBroglie wavelengths being typically much 
smaller than the mean free paths. Of course, the final 
hadron momentum spectra are a product of the time- 
integrated dynamics starting with the initial interpene- 
tration of the two nuclei. The agreement between RQMD 
and experimental data for the final slope parameters im- 
plies that the expansion dynamics in the first ultradense 
stage is modeled reasonably well in this approach. We 
thus believe that the "late" emergence of the baryon flow 
holds true model-independently. 

One might be tempted to look for other explanations of 
the initial "softness" than what is provided by the RQMD 
model. For instance, Hung and Shuryak have put forward 
the idea that the system created in Pb+Pb may be close 
to the so-called softest point of the EOS |3(J. Initial con- 
ditions for hydrodynamical calculations may be tuned to 
get good agreement with the non-strange hadron spectra 
measured in Pb on Pb reactions [ffiflj . However, this hy- 
drodynamical approach fails to explain the mass depen- 
dence of the transverse flow, i.e. the data for S projectiles 
p2[ . On the other side, the RQMD model provides an 
explanation for the systematics from p+A and S+A to 
Pb+Pb Q . According to the RQMD model the trans- 
verse pressure during the early stages is small in all these 
systems, generically a pre-equilibrium signature. 

In summary, we report the results of the analysis of 
the particle transverse momentum distributions from the 
central 158AGeV Pb+Pb collisions. The systematics of 
the baryon spectra indicates that the main component of 
the baryon flow develops only rather late, after most of 
the multi-strange particles have already frozen out. We 
infer from the analysis of RQMD results that character- 
istic energy densities at which the multi-strange hadrons 
freeze out are around 1 GeV/fm 3 for central Pb on Pb col- 
lisions. Presumably, they may be formed at even larger 
densities. 
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